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ABSTRACT 

We present the first clustering results of X-ray selected AGN at z ~ 3. Using Chandra X-ray imaging and 
UVR optical colors from MUSYC photometry in the ECDF-S field, we selected a sample of 58 z 3 AGN 
candidates. From the optical data we also selected 1385 LBG at 2.8 < z < 3.8 with R < 25.5. We performed 
auto-correlation and cross-correlation analyses, and here we present results for the clustering amplitudes and 
dark matter halo masses of each sample. For the LBG we find a correlation length of ro iBC = 6.7 ± 0.5 Mpc, 
implying a bias value of 3.5±0.3 and dark matter (DM) halo masses of log(M„,„,/M0) = 11 .8 ± 0. 1 . The AGN- 
LBG cross-correlation yields roACN-LBC = 8.7 ±1.9 Mpc, implying for AGN at 2.8 < z < 3.8 a bias value of 
5.5±2.0 and DM halo masses of log(M,„,„/MQ) = 12.6;!;y g. Evolution of dark matter halos in the Lambda CDM 
cosmology implies that today these z 3 AGN are found in high mass galaxies with a typical luminosity l'!i2L*. 
Subject headings: galaxies:active,high-redshift — large-scale structure of universe 



1. INTRODUCTION 

There is a wealth of evidence that nuclear supermassive 
black holes play a significant role in the process of galaxy 
formation and evolution. This has become evident in the 
past few years with the discovery of correlations between 
the properties of the massive black holes, th e stellar systems 
that host them, and their dark matter halo s (iMagorrian et al.l 
[T998llGebhardt et al.ll2000HFerraresell2()()2h . 

Luminous quasars have been studied in great detail, with 
systematic spectroscopic studies by the Sloan Digital Sky Sur- 
vey and the 2dF Galaxy Redshift Survey. These objects are 
rare and represent the bright end of the active galactic nuclei 
(AGN) luminosity function. In order to fully understand the 
link between the growth of supermassive black holes (SMBH) 
and galaxy evolution, we need to study young, high-redshift 
galaxies hosting AGN with more typical luminosities. 

One of the basic properties of galaxy populations is 
their clustering strength, but there are few constraints on 
this quantity for fainter, high-redshift AGN. Furthermore, 
there are important disagreements in the literature: e.g. 
galaxy-AGN cross-correlatio n measurements at z ~ 3 by 
lAdelberger & Steidell (l2005bl) (lA05bl) . imply a bias factor 
of 3.9 ± 3.0 for luminous AGN (UV luminosities between - 
30< Mi35o <-25), while recent work by Shen et al. {200% 
(IS07h suggests a value of 9.1 ±0.9 for sources with similar 
luminosities at the same redshift. For galaxies, e.g., measure- 
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ments indic ate a strong luminos ity dependence in the cluster- 
ing length ( Giavalisco & Dickin son 2001 ). On the other hand, 
for A GN there are a handful of claims t hat this is not the case 
(Croo m et al.ll2005l iMvers et al.ll2006l lAdelberger & Steidell 
2005a). Determining the bias with accuracy puts important 
constraints on models of AGN formation and evolution, (e.g., 
iLidz et alJi20()6ilHopkins et al.i 2007) but requires spanning a 
broad range in luminosity and obscuration level. This is chal- 
lenging using optical color selection plus spectroscopy, be- 
cause the low fraction of AGN (3%) found among these can- 
didates (Steidel et al. 2002) demands huge amounts of spec- 
troscopy time . Using X-ray detection plus spectroscopy (e.g., 
ISzokoly et al. 1 12004) provides a more efficient and unbiased 
selection of AGN, since their surface density is higher in X- 
rays than in optical images and obscuration effects are much 
less important. However, the unrestricted redshift range sam- 
pled by this method makes spectroscopic follow-up highly in- 
efficient, and often prevents confirmation of dimmer, high- 
z AGN (R>24). Hence it is difficult to obtain AGN sam- 
ples suitable for clustering studies. There have been sev- 
eral measurements of the spa tial correlation fun c tion of X- 
ray selected AGN at z < 1 dMulhs et al.l j 2004, Gil U et al 
20051. iBasi lakos et al.ll2004llBasilakos e'taLll2005. Yang et aL 
2006, Miv aii et al.l2007l) . with samples ranging from 200-500 
sources. At higher redshifts, the statistics are much poorer, 
and consist of purely optically selected AGN. In this work, we 
constrain the clustering strength of an AGN sample at z 3 
jointly selected by optical and X-ray photometry. We deter- 
mine whether these sources cluster more or less than non- 
active galaxies at this redshift, and discuss their present-day 
descendants. 

We a ssume a ACD M cosmology consistent with WMAP 
results (ISpergel et al.ll200 7) with 0,„ = 0.3, f^A = 0.7, Hq = 70 
km s~' Mpc and erg = 0.8. All quantities are comoving: cor- 
relation lengths scale as h^^, number densities as /i^,, and halo 
masses as /i^g. 



2. OBSERVATIONS 

The MUSYC survey was optimized to study galaxies at z ~ 
3, with imagi ng depths dow n to the spectroscopic limit, U, 
B,V,R^ 26 (iGawiser et al.ll2006b.) . In the ECDF-S field, in 
particular, the 5a limiting (AB) magnitudes achieved in these 
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bandpasses are 26.0, 26.9, 26.4 and 26.4 respectively. The 
main MUSYC catalog is based on the sources detected on the 
combined B VR image, and aperture photometry is perform ed 
on each filter in those positions (see lGawiser et al.ll2006al) . 
The deep Chandra observations of this fie l d have produced 
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four X - ray catalogs: 



(I2003h . iLehmeret al" 



Giacconi et a 
(|2005|), and 



(I2002h. [A lexan der et all 
Virani et a l. (2006). The 



first two comprise 1 Ms of exposures inside the central region 
(CDF-S pr oper), covering an area of «0.1 deg^ (PI R. Giac- 
coni). The [Alexander et al.l catalog has reported flux limits 
of 5.2 X 10"'^ erg cm"^ s"' and 2.8 x 10"'^ erg cm"^ s"' in 
the soft (0.5-2 keV) and hard (2-8 keV) bands, respectively. 
The last two catalogs come from the four «250 ks pointings 
that cover an area of «0.3 deg- around the former field (PI 
N.Brandt). Limiting X-ray fluxes in the extended region are 
1.1 X lO""" erg cm"^ s"' in the soft and 6.7 x lO""" erg cm"^ 
s"' in the hard bands, correspondingly. 

The MUSYC spectroscopic follow-up program carried out 
with Magellan/BaadeH-IMACS has yielded 280 successful 
identifications of z > 2 galaxies. Data were obtained with a 
resolution of R=640 (470 km/s at 5000A) and sHdets of 1.2" 
(P. Lira et al., in prep). The broad wavelength coverage pro- 
vided by IMACS, between 4000 and 9000A, allows the detec- 
tion of Lya and C IV 1549, the two most prominent ultraviolet 
AGN emission lines, at redshifts 2.3 < z < 4.0. 

3. AGN AND LBG SAMPLES 

The full set of unique X-ray counterparts is taken from the 
MUSYC ECDF-S X-ray catalog (Cardamone et al. 2007, 
submitted to ApJ), c onstructed joining the catalog s from 
iGiacconi et aP ( l2002l ). [Alexander et al. (2003), Vira ni et al.1 
(120061) andlL ehmer et al. (2005), using a likelihood procedure 
dBrusa et al.ii200 5) to match sources between the X-ray and 
optical (BVR-selected) catalogs. 

Fig. 1 shows the U -V versus V -R color-color plot for 
sources in the ECDF-S from MUSYC photometry. Lyman 
break color selection ('UV R') corresponds to the region out- 
lined in the upper left side (I Steidel et al.lll996llGawiser et al.l 
|^06a). For this color selection, t/-band fluxes are required to 
be detected at l-cr, and are otherwise set to an upper limit 
equal to their l-cr error This avoids interlopers (typically 
dwarf stars) with uncertain photometry at the cost of incom- 
pleteness in the sample. Additionally, R < 25.5 is required 
to allow for spectroscopic confirmation. Sources presenting a 
drop in the U filter due to intergalactic absorption, along with 
a blue continuum between the V and R filters that rules out 
heavily reddened lower-z objects, are expected to lie at 2.7 
< z < 3.7. For XUVR selection of AGN we require an X-ray 
detection with counterpart UVR colors in the same region and 
we drop the R magnitude limit and the l-cr f/-band require- 
ment, since the extra requirement of X-ray emission already 
rules out most dwarf stars. This procedure yields 1385 LBG 
and 58 AGN candidates. Unobscured AGN and LBG have 
somewhat similar UVR colors, as can be seen in the color- 
color track in Fig. 1, since these colors are mainly determined 
by the intergalactic hydrogen absorption. Furthermore, ob- 
scured AGN are the dominant population among AGN, and 
sin ce the SED of the first is dominated by th eir host galax- 
ies (iTreister et alJl2004lTreister & Urryll2005h . we expect the 
redshift distribution of our AGN sample to differ mildly from 
the LBG distribution. Notice that in the case of significant ob- 
scuration, we don't expect to select AGN hosted by very red 
galaxies, since they will probably not be bright enough in the 
restframe UV. 
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Fig. 1. — U — V versus V — R color-color plot for all the sources brighter 
than R = 25.5 in the MUSYC ECDF-S optical catalog (gray points) and X-ray 
counterparts (open triangles). LBG and ~ 3 AGN selection are inside the 
region marked in solid lines. The color-color tracks of a template AGN and a 
LBG as they are redshifted from z=2.0 to 3.9 are shown in the red solid and 
blue dotted curves, respectively. 
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Fig. 2. — Redshift distribution of all spectroscopically confirmed LBGs 
(white bars) and z 3 AGN (black bars) in MUSYC. The subset of AGN in 
ECDF-S is in gray bars. These distributions have mean and variances: zagn 
= 3.18, a, AGN = 0.39, zlbg = 3.23, ct-^xbg = 0.22. 

From the spectra obtained over the entire 4-field survey 
so far, 131 LBG and 30 AGN were identified at the target 
redshift ~3, and in the ECDF-S in particular, 31 LBG and 
1 1 AGN have been confirmed. While AGN in the ECDF-S 
were directly targeted using a joint X-ray and optical selec- 
tion, in the rest of the MUSYC fields AGN have been dis- 
covered serendipitously among the UVR candidates (see Fig. 
1). Fig. 2 shows the redshift histograms of the confirmed 
objects, which implies X-ray luminosities between 10'*^ and 
10'*^ ergs/s for the entire set of AGN candidates. The Lyman- 
Break (UVR) selection shows a fraction of low-redshift in- 
terloper less than 10%. In the spectra of the XUVR-targeted 
AGN in ECDF-S, we did not observe any interlopers. 

4. CLUSTERING ANALYSIS 

We calculate the clustering of AGN via a two-step method, 
calculating the autocorrelation of LBG first, and then their 
cross-correlation with AGN. The main advantage of this ap- 
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proach is that it improves the statistics significantly, since 
the LBG population is much more numerous and its red- 
shift distribution is similar to the AGN population sampled 
here, the AGN-LBG cross-correlation fu nction less noisy than 
the A GN autocorrelation function (Kauffmann & H aehneltl 
l2002h . Notice that the AGN and LBG that contribute to the 
cross-correlation function will be those that spatially overlap. 
Therefore, our results will reflect the clustering of AGN in the 
redshift range 2.7 < z < 3.8. 

For the UVR selected sources, we calculated llie ang ular 
autocorrelation function using the lLandv & Szalavl(ll993h es- 
timator: 

^ DD{0)-2DR{0)+RR{e) 

""^^^ = me) 

where DD, DR and RR are the data-data, data-random and 
random-random pairs, all normalized to integrate to 1. To 
estimate the angular cross-correlation between UVR-XUVR 
sources, we used the following cross-correlation estimator 
(lCroftetal.ll 19991) : 



DacnDlbc(&) 



1 



(2) 



where now the data-data pairs are between AGN and LBG, 
and data-random between AGN and the optical random cat- 
alog. Notice that this estimator does not use a random cat- 
alog for the X-ray sources. The degradation of the PSF and 
sensitivity as a function of the off-axis angle in the Chandra 
data makes the X-ray survey geometry fairly complicated, but 
these resulting systematic uncertainties can be avoided in this 
cross-correlation scheme. 

The estimates shown above were fitted to (jj{9) = A9~^ - 
lC{A,f3), where the integral constraint factor IC is included 
as part of the model and we assumed the usual power law 
for the true angular correlation function. We kept f3 = 0.8 
fixed, since our S/N does not allow us to put significant con- 
straints on the slope. Fit limits have been set to avoid the 
one-halo term regime on small angular scales (30" for LBG 
and 60" for AGN based on the virial radii of their initially in- 
ferred DM halo masses of 6x lO" and IO'^Mq, respectively) 
and ending at half the field size on the larger scales, in order 
to avoid border effects of the estimator and where the sam- 
pling becomes poor anyway. The fit to the LBG autocorre- 
lation function gave Albc = 2.2 ±0.3 arcsec^ and the fit to 
the cross-correlation function Aacn-lbc = 2.9 ±1.1 arcsec^, 
with reduced values of 1.24 and 1.21, respectively. The er- 
rors in these parameters were calculated using Ax^ and corre- 
spond to 1-(T confidence level with 1 parameter. Fig. 3 shows 
the measured angular correlation functions with their corre- 
sponding best fit models. Binning Lo{0)/a^{9) over the entire 
fitting range gives estimates for the total signal-to-noise ra- 
tios of these measurements of 10.2 and 3.9, correspondingly. 
From montecarlo realizations, we estimated that the probabil- 
ities of obtaining clustering this high from unclustered popu- 
lations of the same sizes are <0. 1 % for the LBG autocorrela- 
tion and 1 % for the AGN-LBG cross-correlation. 

Using smooth fits to the redshift distributions of the con- 
firmed sources shown in Fig. 2, we deprojected the values 
obtained for the angular LBG autocorrelation and LBG- AGN 
cross-correlation amplitude using Limber's equation. We ob- 
tained spatial correlation lengths of roxBC = 6.7 ±0.5 Mpc and 
rQ AGN-LBG = 8.7 ± 1 .9 Mpc respectively, where the correlation 
length is defined by the usual fitting form ^(r) = {r/ro)~^ for 
the spatial correlation function. The random errors introduced 
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Fig. 3. — Top: Angulai' autocon'elation function for UVR selected sources, 
i.e. Lyman-break galaxies at z ~ 3. The solid line is the best fit to the data, 
^lbg{S) = 2.2 ± 0.3 9"" **. The point marked with a star is not included in the 
fit and is shown only for reference. Bottom: Analogous plot for the angular 
cross-correlation function with XUVR sources, i.e., AGN at z ~ 3. The best 
fit corresponds to tjJagn-lbg{6) = 2.9 ± 1 . 1 fi""'* . 



TABLE 1 
Summary of DM halo properties 





bias 


log(M„„„) 


nimlo (Mpc ^) 


LBG 


3.5 ±0.3 


11.8±0.15 


9±4x K)-^ 


AGN 


5.5 ±2 


12 6+"-' 


10"' - 10"* 



relations b]^g = ajj^^ 



where a^x is the variance in spheres of radius 8/1 



Note. — Confidence intervals/ranges are l-cr. Masses are given in loga- 
rithm of solar masses. 

in the deprojection are included in the error budget by approx- 
imating the redshift distributions as Gaussians and propagat- 
ing the standard errors in their means and variances. 

To calculate the bias factor of the AGN sample in the ap- 
proximation of linear, constant bias, we used the elementary 

agnblbg = (^iMgn-lbgl ^%.dnv 

Mpc 

in the random field described by the corresponding auto- or 
cross-correlation function. Although A//? and ro/7 are highly 
degenerate parameter pairs, the bias factors are robust to vari- 
ations in the slope (3 (or 7). Therefore, they are our pre- 
ferred quantity for comparison with the literature. The bias, 
number densities and masses shown in Table 1 are calculated 
using the ellipsoidal co Uapse model extension o f the Press- 
Schechter form alism bylSheth & TormenI (1 19991) and the ex- 
tension of the iMo & White! {\99(h formalism by I Sheth et all 
(2001). These quantities were evaluated at the mean redshifts 
of the confirmed samples, both very close to z=3.2. 



5. DISCUSSION 

We have measured the autocorrelation strength of z=3 LBG 
and the cross-correlation between these galaxies and AGN 
selected using both X-ray and optical data. For the LBG 
sample, we obtained a bias factor of 3.5±0 .3. This is some- 
what h igher th an the value 2.8±0 . 3 foun d bv lAdelberger et al] 
(l2005h at z = 3. lHildebrandt etalJ (l2007h obtained a bias value 
of 3.2±0.2 for an equivalent LBG population, consistent with 
ours. 

From the AGN-LBG cross-correlation and the bias calcu- 
lated for the LBG, we have deduced a bias factor of 5. 5 ±2 for 
our AGN sample. The active galaxies targeted in this study 
appear to cluster more than star-forming galaxies with similar 
restframe-UV colors. This is consistent with cosmic downsiz- 
ing of AGN, implying that typical SMBH tend to sit in more 
massive galaxies than the "normal" galaxy population. We 
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Fig. 4 . — Comparison of AGN bias factors at 7, ' 
(square), §07 (triangle) and lAOSbl (stars). 
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need greater statistics to confirm this result, since the values 
are consistent within the uncertainties. 

A05b performed this same calculation in an optically se- 
lected and spectroscopically confirmed sample of 79 AGN be- 
tween 1.6 < z < 3.7, dividing the AGN by UV luminosity into 
bright (25 sources with -30< M1350 <-25) a nd dim (54 sources 
with -25< M1350 <-19) samples. Although 'A05b' does not re- 
port the autocorrelation length of the galaxy sample used to 
compute these values, we approximate it by avera ging the re- 
sult fo r LBG and BX galaxies presented in Adelberger et al.l 
(l2005h . obtaining a bias factor of 2.6±0.3. F rom th is we infer 
bias factors of 3.9±3.0 and 4.7±1.7 for the 'A05b' bright and 
dim AGN samples. Our AGN sample has an UV magnitude 
range between -26 and -20, almost identical to the A05b faint 
AGN set, and showing the same clustering strength. In Fig. 4 
these result s are compared to the present estimate and to that 
obtained by'SOT'for ^2250 SDSS quasars at 2.9< z <3.5. 

To estimate the dark matter halo mass of the typical de- 
scendant of the halos that host these z galaxy sets at the 



present time, we calculated the mode and width of the condi- 
tional probability di stribution of the expected mass z = (see 
iHamana et al.l|20 d^ and references therein). We took the me- 
dian mass of the halo populations at z = 3.2 as the representa- 
tive value, and we report the bias values corresponding to the 
maximum likelihood z=0 progenitor, with \ -a uncertainties 
corresponding to halo masses whose likelihood is reduced by 
a factor exp(-\ /I) = 0.607. For the LBG bias of 3.5 at z = 3.2, 
we obtain a bias of 1.3;!;o j at z = 0. In the nearby universe 
dZehavi et al.l l2005l) . this corresponds to the clustering of a 
somewhat massive galaxy, with L = . On the other 

hand, our AGN sample, with a bias of 5.5 at z = 3.2 would 
have a typical ha lo with a bias factor of 2.0^o 3 at present. Ex- 
trapolating from lZehavi et alJ (l2005h . this corresponds to the 
clustering of galaxies having It^L*, which at the present time 
are most likely located in groups an d galax y clusters. 

This confirms the result found by iA05bl for faint, optically 
selected AGN. Improving our understanding of the AGN- 
galaxy connectio n at z 3 and con straining models such as 
iLidz et alJ 42006") and lHopkins et all (12007) for AGN cluster- 
ing, requires resolving the discrepancy in the bias estimates 
for bright AGN seen in Fig. 4 improving the statistics. We 
have show that targeting high-redshift AGN for clustering 
analyses can be done very efficiently by means of deep optical 
and X-ray imaging, and for that reason, future surveys with 
deep X-ray coverage will be ideal for obtaining large sam- 
ples of active galaxies in restricted redshift ranges, suitable 
for clustering studies. Our method can be applied to these 
surveys to obtain enough X-ray-selected AGN to reduce the 
current uncertainties. 

H. F was supported by MECESUP project UCH0118, 
Andes Fundation fellowship C- 13798, ALMA fellowship 
31060003 and Fondecyt project 1040719. 
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